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Near a parity breaking front bifurcation, small perturbations may reverse the propagation direc¬ 
tion of fronts. Often this results in nonsteady asymptotic motion such as breathing and domain 
breakup. Exploiting the time scale differences of an activator-inhibitor model and the proximity to 
the front bifurcation, we derive equations of motion for planar and circular fronts. The equations 
involve a translational degree of freedom and an order parameter describing transitions between left 
and right propagating fronts. Perturbations, such as a space dependent advective held or uniform 
curvature (axisymmetric spots), couple these two degrees of freedom. In both cases this leads to a 
transition from stationary to oscillating fronts as the parity breaking bifurcation is approached. For 
axisymmetric spots, two additional dynamic behaviors are found: rebound and collapse. 


I. INTRODUCTION 

Pattern dynamics in reaction-diffusion systems often 
involve nonsteady front motions. These motions can be 
driven by curvature , front interactions [pH], con¬ 
vective instabilities l5|% and external fields In 

some cases fronts reverse their direction of propagation, 
as for example in breathing pulses , where the re¬ 

versal is periodic in time, and nucleation of spiral-vortex 
pairs^here the reversal is local along the extended front 
line [ll|-||. 

Earlier studies demonstrated that front reversals, as 
described above, become feasible near a nonequilibrium 
Ising-Bloch (NIB) bifurcation that is, a parity 

breaking bifurcation where a single stationary front loses 
stability to a pair of new, counterpropagating fronts. The 
reversal phenomenon can be regarded as a dynamic tran¬ 
sition between the left and right propagating fronts that 
appear beyond the front bifurcation. It is induced by 
intrinsic perturbations, like curvature and front interac¬ 
tions , or weak space dependent external fields . 

Since the left and right propagating fronts differ in inter¬ 
nal structure [^| such a transition involves a new degree 
of freedom in addition to the translation mode: the or¬ 
der parameter associated with the bifurcation |^]. The 
effect of the perturbations is to couple these two degrees 
of freedom in a way that allows for front reversal. 

Our objective is to derive equations for front motion 
that capture front reversal. Progress toward that goal 
has already been made in Ref. 0,^ for a nondiffus¬ 


ing inhibitor. Since inhibitor diffusion is essential for 
spontaneous front reversals induced by curvature, in this 
work we study the more difficult case of a diffusing in¬ 
hibitor. This calls for a different approach as described 
in Sections HI and IV. In Section ^ we study front re¬ 
versals induced by two types of perturbations of planar 
fronts, an external advective field and uniform curvature. 
In this case only planar and circular fronts are studied; 
the derivation of the more general equations for nonuni- 
formly curved fronts can be found in Ref. [0. Some of 
the results presented here have been briefly reported in 


II. REACTION-DIFFUSION MODEL 

The model we consider is an activator-inhibitor 
reaction-diffusion system describing a bistable medium. 
Models of this type have been studied in various physi¬ 
cal and chemical contexts [^]j2|0-0|. The specific form 
chosen here is 

ut = e~^{u - - v)5~^Uoi;x , (1) 

vt = u - aiv Vxx ■ 

The variables u and v are scalar real fields repre¬ 
senting the activator and the inhibitor, respectively, 
with the subscripts x and t denoting partial derivatives 
with respect to these variables. For ai > 1 the sys¬ 
tem (|i|) has two stationary uniform states, (m±,u±) = 
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FIG. 1. The NIB (or nonequilibrium Ising-BIoch) bifurca¬ 
tion and internal structure of front solutions. The pitchfork 
diagram represents the speed of front solutions vs the param¬ 
eter rj. For rj > rjc, the Ising front is the single solution and 
Vf, the order parameter representing the value of the v field 
at the front position u = 0, is zero. Beyond the bifurcation, 
r] < rjc, a. pair of counterpropagating Bloch fronts appears. 
The order parameter is negative (positive) for rightward 
(leftward) propagating fronts. 

^±y/l — 1/ai, — 1/ai^ . Note that the parity 

symmetry {u,v) —>• (—u,—v) of Eqs. ([^) is reflected in 
these solutions. Generalizations of Eqs. ([|) to nonsym- 
metric forms in one and two space dimensions will be 
considered in Section ^ 

In addition to the spatially uniform solutions there are 
also front solutions connecting them. In the following 
we will consider front solutions that connect (u+,u+) at 
X = — oo to (u_,u_) at a; = oo. The number and type 
of these front solutions is determined by the two param¬ 
eters e and 6. For rj > rjc = Sj2'/2q^, where rj = \/eS 
and = ai + 1/2, a single stable stationary (Ising) front 
solution exists. This solution loses stability in a pitchfork 
bifurcation, at 77 = 77 c, to a pair of counterprop^ating 
(Bloch) front solutions as shown in Fig. |l[ The 

two Bloch front solutions differ not only in their prop¬ 
agation direction but in their internal structure. This 
difference can be represented by an order parameter as¬ 
sociated with the bifurcation, which for /i = e/(5 <C 1 can 
be taken to be the value u/ of the v field at the front 
position. For simplicity we define the front position to 
be at M = 0. With this choice, Vf=0 for the Ising front 
(the inset for 77 > 77 c in Fig. p. Beyond the front bifur¬ 
cation, 7 ;/ is nonzero and the sign indicates the direction 
of front propagation: Vf < 0 for the front propagating to 
the right and Vf > 0 for the front propagating to the left 
(77 < 77 c in Fig. p. 


III. FORMULATION OF THE FREE BOUNDARY 
PROBLEM 

In the following we confine ourselves to the region 
e <C 1 , i5 cx e~^ and we choose d values such that 
ed ~ 0{ril). The small parameter e allows the use of 
singular perturbation methods to study front solutions 
to Eqs. (|^). The first step is to transform to a moving 
coordinate frame, x ^ r = x — Xf(t), t ^ t where Xf is 
the position of the front. In this frame Eqs. ( 0 ) become 

Ut — if Ur = e~^{u — — v) + S~^Urr , (2) 

Vt — XfVr = U — aiV + Vrr , 

where the dot over x / denotes the derivative with respect 
to t. The front solution, u{r,t), v{r,t), is characterized 
by a strong variation of the u field over a distance of 
order = \fe[6. Stretching the spatial coordinate, 
z = r/y/JI, to expand this region Eqs. (||) become 

e{ut - ZfUz) = u-- V+ Uzz , (3) 

fL (vt - ZfVz - U + aiv) = Vzz , 

where Zf = Xf f y/JI and we recall that /i oc e^. Expanding 
in e 

u = uq + eui + e^U 2 + ..., 

V = Vo + f-Vi + €^V2 + ..., 

and inserting into we find at order unity the front 
solution 

Uq = — tanh(z/-\/ 2 ), uq = 0 . 

Collecting terms of order e gives 

£ui =vi- ZfUoz , £ = + 1 - 3uo , (4) 

where vi is a yet undetermined function of time. Since 
i^UQz — 0 ; 

solvability of (||) requires 

The narrow front region becomes infinitely thin in the 
limit e —> 0. Therefore, v{t) may be associated with the 
value of v{r,t) at the front position, that is v{Q,t). With 
this notation the leading order relation is 

*/=(5) 

Away from Xf, u — — v ^ 0(e), and u varies on the 

same time and length scales as v. Going back to Eqs. (^) 
we find at order unity 
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Vt — XfVr = U+{v) — aiV + Vrr , r <0 , ( 6 ) 

Vt — XfVr = U-{v) — aiV + Vrr , X > 0 , 

where u± (v) are the outer solution branches of the cubic 
equation u — u^ — v = 0. For oi sufficiently large we may 
linearize the branches u±(y) around v = 0 

« ±1 — u/2 . (7) 


Note that = —< 7 “^ exp(—q|r|) is not differentiable at 
r = 0 . 

In terms of the deviation from the stationary solution, 
V = V — , Eqs. are 

Vt + q^V - Vrr = - ^v{0,t){vr + ) , (12) 

?7V2 

u(±oo) = 0, 


Substituting the linearization (|^) and the relation from 
the inner problem (||) into produces the free boundary 
problem 


r < 0, 


r > 0, 
( 8 ) 
(9) 

where the square brackets in denote jumps across the 
free boundary. The solution to (||) leads to a dynamic 
equation for Vf{t) = u(0,f), the value of the inhibitor at 
the front position Xf{t), which will complement Eq. 


Vt+q^V—Vrr = 1- j=v{0,t)l 

77V2 

v{—oo,t) = 


Vt+q^V — Vrr = — 1 - -j=v{Q,t)l 

77V2 

v{oo,t) = v-~—q~‘^ 


b]r=0 “ ^ , 


IV. SOLUTION OF THE FREE BOUNDARY 
PROBLEM 


Near the front bifurcation (77 close to rjc), the front 
speed, c, is small and propagating front solutions can 
be expanded as power series in c around the stationary 
front solution. The stationary front solution satisfies the 
boundary value problem 


Vrr — q^V + 1=0 

v{—oo) = q~^ 

u(0) = 0 


r < 0, 


c\ 

Vrr — q V —1 = 0 

v{oo) = —q~^ * 

u(0) = 0 


r > 0 , 


with = 0. The solution to (p^ is 


v^°\r) =q ^(1 —e'?’'), ?" < 0 , 

— 1) , r > 0 . 


( 10 ) 


( 11 ) 


We seek propagating solutions of ( 0 ) that involve two 
time scales, the original time t and a slow time T = c^t. 
The slow time dependence is a result of the slow evolution 
of u/ (the value of the inhibitor at the front position) near 
the front bifurcation. It is easy to show that to linear 
order Vf oc {rjc — 'q)vf and for a pitchfork bifurcation 
r]c — 1 ] (X c^, hence the scale. Expanding v(r,t,T) in 
powers of c and rj in powers of (expecting a pitchfork 
bifurcation), 


;(r,t,r) = ^c%(")(r,t,T), 


(13) 


n—1 


rj = r]c- c^rji + c^r ]2 + ..., 


(14) 


and inserting in (p^) gives the sequence of equations 




where 


+ = _p("), ri = l,2,3, (15) 






„(o) 




V2r]< 

^ 1^(1) ^(1) + ^,(2) ^(0) 




qr=0‘' 


„(3) _ ^,(1) , ^,(1) ^(0) 


(16a) 

(16b) 

(16c) 


V2r]c 


(1) i;(2) . ^(2) Jl) , ^(3) +0) 


qr=o 


|r= 0 ‘- 


Equation (|l^ can be solved using an appropriate 
Green’s function and assuming that the front motion is 
independent of the fast time scale t as t ^ 00 |2^ . A sim¬ 
plified derivation of the solution follows from the gradient 
nature of (p^ when the source term p^'^\r,T) becomes 


independent of t. Eor then v. 


(n) 


0 as f > 00 for any r 


and we can look for stationary {t independent) solutions. 
Consider first Inserting (1^) in Vrr'^ — q^v^^^ = 

and solving for we obtain 


.(1) 


(r,T) = 


2\/2r]c 


,uW(0,r)F(r), 


(17) 


where 


F{r) = {1 — qr)e'^^ , r<0, 

F{r) = (1 + qr)e~‘^'~ , r > 0 . 
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Setting r = 0 we find 


Vc = 


2V2q^ ■ 


(18) 


The critical value rjc = rj{c = 0) determines the bifurca¬ 
tion point where the propagating Bloch front solutions 
coincide with the stationary Ising front. The expres¬ 
sion (|^ is the same as the one derived earlier using 
a different method 15 
Using (0) to solve Vrr — 

i;(2) (r, T) = ( 0 , T) + (0, F{r), (19) 


m 


QVc 
6^2 ■ 


( 22 ) 


Equations ( |2^ ) and (|^) provide the leading order form 
of the front bifurcation diagram 


c 


2 



v), 


which coincides for small c with the earlier result [6,9], 
= 4g^(r7^ — p^) /r?^. 

Multiplying ( pi|) by c and using the expansions ([l^ ) 
and (|^) gives the equation of motion for propagating 
fronts 


and using both ( 0 ) and ( p^ ) in viV — = p(3) we 

obtain 


■p){r,T) = z;U)(0,T)e9’'-h A+re«’' 
- - C'+r^e'^’', 


r < 0 , 


= nU)(0,T)e-9''-h A.re-'?’' 


where 


r > 0 , 


( 20 ) 


A± = g3^U)(o,T)wU)(0,T) 

± l4i)(0,T) + igV')(0,r)3 

- —uU)(0,r)-g?;U)(0,T) 

Vc 

B± = iuW (0, T) ± g4^;U) (0, T)uU) (0, T), 

C±=±igVi)(0,T)U 

6 


3 

= S<’’“ “ ■ 3tf”' ’ 

, (23b) 

where the slow time derivative oi Vf is expressed in terms 
of a fast time derivative (vf = c^Vfj,). 

According to Eqs. (^ the dynamics of a propagat¬ 
ing front involve two degrees of freedom: a translational 
degree of freedom, xj, determining the front position, 
and an order parameter, n/, determining the direction of 
propagation. The latter has not been appreciated enough 
since most works to date ||,^ have focused on conditions 
far from the front bifurcation. In that case the two sta- 
tionary states, u/ = ± [A'/^{rjc — ??)/3g] , represent¬ 

ing fronts propagating in opposite directions, are highly 
stable. Close to the front bifurcation, however, the eigen¬ 
value associated with these states, A = —2y/2{r]c — ri)/q, 
approaches zero and small disturbances can drive the sys¬ 
tem from one state to another, thereby inducing front 
reversals. 



V. FRONT REVERSAL: 
OSCILLATIONS AND REBOUND 


Application of the (no) jump condition = 0 

leads to 

(0, T) = (0, T) - (0, Tf . (21) 

Equation ( 0 ) still contains an unspecified parameter, 
rji- To identify 771 recall that c is the speed of a front 
propagating at constant velocity. From (||) 

\if \ = ■y|^c|u(i)(0,T)| -hC>(c2). 

Identifying ji/j with c gives = 2t]1/9 for a front 

propagating at constant speed. This value of ?;U) should 
coincide with the nontrivial stationary solution of (|^), 
z;U) = 4-\/27?i/3g. Comparing the two expressions gives 


Equations (^ describe the motion of a freely prop¬ 
agating front in a uniform medium. In this Section we 
show how two different perturbations affect front prop¬ 
agation. The first is the addition of a space dependent 
advective field to the v equation in the original system 
(P. This type of differential advection appears for ex¬ 
ample in chemical reactions involving ionic species that 
are subjected to electric fields |^,Q. The second is the 
intrinsic perturbation of uniform curvature variations on 
the propagation of two dimensional fronts. Both pertur¬ 
bations lead to a coupling of the two degrees of freedom 
in the order parameter Eqs. ( p^ and allow for the non¬ 
steady asymptotic motion of fronts. 
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A. Space dependent advective field 


To study the effect of an external advective field J we 
add the term Jvx to the v equation in (|^), 


-i„ 


Ut = e ^{u — — v) + S 

vt = u- aiv - tto + Jvx + Vx 


(24) 


The small parameter ag is also introduced to break the 
parity symmetry of ( 0 )- For simplicity we consider a lin¬ 
ear spatial profile, J = —ax, 0 < a <C 1. Proceeding 
as before, the inne r reg ion analysis remains unchanged 
and culminates in ( 23b ). The outer analysis leads to the 
additional terms, —a(r + Xf)vr — ag, on the right hand 
side of both partial differential equations in (||) . Assum¬ 
ing a = agc^ and ag = aggc^ , wh ere ao and agg are 
of or der unity, Eqs. (1^) and ( |l6b| ) remain unchanged, 
but (16c) acquires on the right hand side two additional 
terms: ao(r-|-cc f)u r°^ +aoo- As a result the order param¬ 
eter equations (E^ are changed to 


V2. . 3^2 


;ag . 


Xf = 


vV2 


Vf 


(25a) 

(25b) 



FIG. 2. Front position, Xf, vs time for an oscillating front. 
The solid line represents the solution to the order parameter 
equations (^^ and the diamonds are from the numerical solu¬ 
tion of the original partial differential equations (^). Param¬ 
eters: ai = 3.0, e = 0.01, S = 2.77, ao = 0, a = 0.005. 


Notice that the introduction of a space dependent ad¬ 
vective field couples the two degrees of freedom, Vf and 
Xf. This coupling affects the front behavior in two sig¬ 
nificant ways: for p > (and ag ^ 0) it stabilizes a 
propagating front at a fixed position, Xf = 2qagla, and 
for r; < ? 7 c it induces oscillations between the counter- 
propagating fronts. The frequency of oscillations close to 
the Hopf bifurcation at rj = rjc is 

w=-^(7Va. (26) 

We tested the validity of Eqs. ( ^^ by numerically inte¬ 
grating the original system (j^) and comparing the oscil¬ 
lating front solutions with those of (p^. The agreement 
as Fig. H shows is very good. In Fig. || we plotted the 
frequency of front oscillations vs the field gradient ac¬ 
cording to d^ ) and as obtained from (|I). Again, the 
agreement is excellent, and remains good even for c of 
order unity. Note that in the inner analysis we neglected 
contributions of G(e^) to v at the front position, while 
in the outer analysis we kept terms to G(c^). A quanti¬ 
tative comparison as described above therefore requires 
that c is much larger than 

B. Uniform curvature 

In two space dimensions the reaction-diffusion sys¬ 
tem (|^ becomes 



FIG. 3. A log-log plot of the oscillation frequency, uo, vs 
the external field gradient, a. The solid line is the relation of 
Eq. (M) and the diamonds represent numerical solutions of 
Eqs. ®). Parameters: ai = 3.0, e = 0.01, S = 2.77, ao = 0. 
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Ut = e ^{u — u^ — v) + 6 (27) 

Vt = u — aiv — ao + , 


where the small parameter oq has been added again to 
break the parity symmetry of #■ In addition to planar 
front solutions there are now new types of solutions in¬ 
cluding fronts with uniform curvature (circular fronts or 
spots). These spots may be stationary or, for parameters 
near a NIB bifurcation, may collapse, expand indefinitely, 
or oscillate periodically in time. 

To derive equations for the motion of circular fronts, 
the first step is to transform into polar coordinates, r = 
p — Pf(t), that move with the front. In this frame and 
assuming the radius of curvature Pf is much larger than 
the front width, Eqs. ( |^ ) are 

Ut — {pf + S~^K)Ur = €~^{u — — v) + S~^Urr , 

Vt — {pf + K)Vr = u — aiV — Oq -|- Vrr , 

where k = pj^ is the front curvature. As before we 
assume e <C 1 and 5 oc and use singular perturbation 
theory. 

Analysis of the inner, or front, region yields a relation 
analogous to ( |23b| ) 

pj: + S~^K = --^Vf. (28) 

In the outer region, instead of (||) we must solve 

3 

Vt + q^V — Vrr = ±1-p1’(0, t)Vr + P , 

■qy/2 


where P = (1 — S~^)vr/pf — ag. Assuming that P is a 
small perturbation (of order |cp) and proceeding as in 
Section 0 we obtain the order parameter equation 


Vf = —[lie 

mi 



2 ( 1 -^-^) 
3q Pf 



(29) 


Writing Eqs. ( p^ and ( p^ terms of the curvature k = 
pj^ gives the equations 


3 2 4 

-/ = - ^(1 - 

k= ^ VfK^ + S~^K^ , (30b) 

that describe the dynamics of large circular spots. The 
introduction of curvature couples the two equations. 
Equations ( p^ for planar fronts are decoupled and only 
describe the relaxation to steadily propagating fronts. 
The equations for circular fronts additionally allow front 
reversals and nonsteady asymptotic motion like oscilla¬ 
tions. 

Consider first the fixed point solutions obtained by 
the intersections of the linear nullclines k = 0 and k = 


— [3d/r]\/2)vf of (^) with the cubic nullcline of (30a). 
The solutions corresponding to k = 0 describe planar 
fronts propagating at constant velocities. Solutions with 
positive and negative Vf values pertain to down states 
invading up states and up states invading down states, 
respectively. The number of k = 0 solutions varies with 
77 . Below the front bifurcation, [i] > ? 7 c(ao)], there is a sin¬ 
gle intersection point representing an Ising front as shown 
by the thin lines in Figs. ^ and ^b. Beyond the front 
bifurcation, [rj < r]c[o-o)], two more intersection points 
appear corresponding a stable and unstable pair of pla¬ 
nar front solutions (Fig. |c). The fixed point solutions 
for K 0 represent a circular fronts. For oq < 0 they de¬ 
scribe spots of up an state domain and for oq > 0 spots 
of a down state domain. For d > 1, depending on the 
choice of e, these fixed points may or may not be stable. 
For 5 < 1, all the k 7 ^ 0 fixed points are unstable. 

Figure m shows three different possibilities for the dy¬ 
namics of circular fronts. The thick trajectories represent 
dynamics computed by numerical solution of the coupled 
equations d^). The initial conditions correspond to a 
large shrinking up state spot. Far into the Ising regime 
(Fig. |a) the initial spot converges to a stationary spot. 
Moving closer to the front bifurcation and past a critical 
77 value, PH > ? 7 c(ao)j a Hopf bifurcation to a breath¬ 
ing spot occurs (Fig. ^). Crossing the front bifurcation, 
77 < 77 c(ao), the spot rebounds, i.e. the shrinking spot 
reaches a minimal size and expands again indefinitely 
(Fig. |c). For larger |ao|, there is another possibility for 
the dynamics of spots. In this case, shown in Fig. ||, 
the amplitude of oscillations grows in time until the spot 
eventually collapses as the curvature diverges to infinity. 

All the three behaviors discussed above have been ob¬ 
served in direct numerical solutions of (|^). The quan¬ 
titative accuracy of the order parameter equations was 
tested by computing numerical solutions to the circularly 
symmetric version of Eqs. (H) 


Ut = e ^{u — U^ — v)-\ - Ur + S ^Urr , (31) 

r 

1 

Vt = u — aiV — Oo H- Vr + Vrr ■ 

r 

and comparing them to solutions to Eqs. ( ^ ) for spot 
dynamics. Spot solutions of ( |^ produce the same qual¬ 
itative behavior as the pair of coupled equations for the 
spot dynamics. When the parameters are chosen to sat¬ 
isfy the assumptions made in the derivation of (|^), there 
is also quantitative agreement between the two solutions. 
Figure || shows the curvature of an oscillatiM spot as a 
function of time computed using both Eqs. (|^ ) and (^). 
The two solutions agree within an accuracy of approxi¬ 
mately 1% for the amplitude and 2% for the phase. 

In addition to the oscillatory instability spot so¬ 
lutions may also be unstable to transverse perturba¬ 
tions [p5|-p^. Numerical solutions of the fully two- 
dimensional model (|2^) show that for the parameters of 
Fig. 1^ spots are unstable and form nonuniformly curved 
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FIG. 4. Three types of solutions to the order parameter 
equations (^) starting with initial conditions representing a 
large shrinking spot. The thin lines are the isoclines and the 
thick lines are numerically computed trajectories, (a) Con¬ 
vergence to a stationary spot (e = .0063). (b) An oscillating 
spot (e = .006). (c) Spot rebound and expansion of the spot 
to infinite size (e = .0052). Parameters: ai = 4.0, ag = —0.01, 
and S = 2.0. 



FIG. 5. A trajectory of the order parameter equations 
for a spot that oscillates with growing amplitude until collapse 
(the curvature k diverges to infinity). Parameters: ai = 4.0, 
ag = —0.1, e = 0.006, and S = 2.0. 



FIG. 6. An oscillating circular spot solution. The 
solid line is the solution of the order parameter equa¬ 
tions (|30|), and the diamonds represent the spot curva¬ 
ture vs time from the numerical solution to the circularly 
symmetric equations (^). The equation parameters are 
£ = 0.006, 5 = 2.0, ai = 4.0, ag = —0.01. 
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fronts leading to a labyrinthine pattern. Since the order 
parameter equations derived here apply only for the case 
when the spots do not break perfect circular symmetry, 
for this choice of parameters they only capture the dy¬ 
namics of the circular spot during the initial evolution. 
Order parameter equations for the dynamics of nonuni- 
formly curved fronts are presented in Ref. |2^. 


reports on the observation of stationary, breathing, and 
rebounding spots. Interaction between fronts may sim¬ 
ilarly lead to stationary, oscillating, and collapsing do¬ 
mains MM- Recent experiments on the Ferrocyanide- 
lodate-Sulfite reaction show small oscillating chemical 
spots away from the reactor boundary that are most 
likely due to front interactions and/or curvature 


VI. CONCLUSION 


ACKNOWLEDGMENTS 


We derived the equations that govern the dynamics of 
planar fronts in bistable systems near a parity breaking 
front bifurcation (the NIB bifurcation). In this case the 
context is an activator-inhibitor model, but the normal 
form 


X = C, 

C = {ac- a)C - /3C^ 


(32) 


is general. Here, X is the front position, C is the front 
velocity, and is a critical parameter value for which 
a NIB bifurcation occurs. Similar equations should ap¬ 
ply, for example, to liquid crystals subjected to rotat¬ 
ing magnetic fields p8|-^ and have also been proposed 
in the context of parity breaking traveling-wave bifurca¬ 
tions 

Uniform front curvature, or space dependent external 
fields, couple the two degrees of freedom, X and C, and 
allow nonsteady asymptotic behavior. The coupled sets 
of equations (|^) and ( [30| ) exhibit Hopf bifurcations from 
stationary to oscillating fronts (breathing spots). Equa¬ 
tions ( p0| ) exhibit two additional behaviors pertaining to 
rebounding and collapsing spots in the full equations. 

Curvature effects on front dynamics near a NIB bifur¬ 
cation were also studied in Refs. using a “qua¬ 

sistatic” approximation. This approximation, where the 
front velocity is assumed to adiabatically follow slow cur¬ 
vature variations yields an algebraic relation be¬ 

tween the front velocity C and its curvature k. As the 
bifurcation is approached the C—k relation becomes mul¬ 
tivalued, or hysteretic. The multivalued relations cor¬ 
rectly predict spontaneous front transitions induced by 
curvature |^,^, but cannot describe dynamics during 
front transitions. Differential order parameter equations, 
like (|^), give a more accurate characterization of the 
dynamics. These differential equations reduce to an al¬ 
gebraic C — K relation when the time scale of front transi¬ 
tions becomes much shorter than the time scale of curva¬ 
ture changes. Such a condition is realized, for example, 
with very large spots away from the front bifurcation. 
Then the right hand side of (30a) can be set to zero, an 
expression which together with ( |28| ) gives an algebraic 
C — K relation where C = pf. 

The phenomena of breathing, oscillating, and collaps¬ 
ing spots appear to be quite general and can be induced 
by other perturbations that couple X and C. Ref. Q, 
which studies the effect of boundaries on spot dynamics. 
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